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OVERVIEW 


The program objectives are fully defined in the original proposal entitled “Program of 
Research in Flight Dynamics in GW at NASA Langley Research Center,” which was originated 
March 20, 1975, and in the renewals of the research program from January 1, 2003 to September 
30,2005. 

The program in its present form includes three major topics: 

1 . the improvement of existing methods and development of new methods for wind 
tunnel and flight data analysis, 

2. the application of these methods to wind tunnel and flight test data obtained from 
advanced airplanes, 

3. the correlation of flight results with wind tunnel measurements, and theoretical 
predictions. 

The Principal Investigator of the program is Dr. Vladislav Klein. Two Graduate 

Research Scholar Assistants (M. M. Eissa and N. M. Szyba) also participated in the 

program. 



SPECIFIC DEVELOPMENTS 


Specific developments in Aerospace during the period January 1 , 2003 through 
September 30, 2005: 

1. Data Analysis of highly swept delta wing aircraft from wind and water tunnel data 

a) Wind tunnel experiments were performed on the 1 8% scale model of the 
F-16XL aircraft in the NASA Langley Research Center 14x22 Foot 
Subsonic Wind Tunnel in 1997 and 2004 and included static tests and 
dynamic forced oscillation tests in roll, pitch and yaw angles. Data were 
collected with varying mounting techniques for both static and dynamic 
tests. Data from static tests were used to compute non-dimensional 
aerodynamic derivatives with respect to angle of attack and sideslip angle. 
Results varied for different mounts and Reynolds numbers. 

The measured data from forced oscillation tests were obtained at different 
angles of attack, amplitudes, Reynolds numbers, and frequencies. Some 
of the tests were performed with initial offsets in either roll or yaw angle, 
for roll or yaw oscillations, respectively. Harmonic analysis was applied 
for obtaining the in-phase and out-of-phase components of aerodynamic 
coefficients. Responses varied with amplitude, Reynolds number, 
frequency, and mounting method. In general, decreased model accuracy 
was observed at angles of attack greater than 30°. Increased higher 
harmonic in the math model improved accuracy for large amplitude and 
offset data. Harmonic analysis as well as power spectral analysis 
indicated an increased measurement noise at these low amplitude 
oscillations. 

Frequency dependence was observed mostly between 20° and 50° for all 
components in pitch, roll and yaw oscillations. For frequency dependent 
data a two-step linear regression algorithm was used to obtained 
parameters in an unsteady aerodynamic model. The static coefficients 
computed from forced oscillation tests were compared to static 
coefficients measured from static tests. Axial and normal force 
coefficients from pitch oscillations had high model accuracy. Rolling 
moments coefficient had the highest accuracy for roll oscillations when 
the linear model was used. There was some variation in results between 
small and large amplitudes. 30° amplitude oscillations at zero offset had 
high model accuracy, but the necessity for higher order terms in harmonic 
analysis brought into question the applicability of the relationships 
between the first order Fourier coefficients and the aerodynamic 
coefficients. The variation or the non-dimensional time constant showed 
an increased unsteady effect between 35° and 45° angle of attack. Very 
few yaw oscillations had acceptable unsteady model accuracy. 


The results from wind tunnel data are summarized in NASA TM and M.S. 
Thesis (see Publications 1 and 2). Some examples are included in Figures 
1 to 4. 

b) The problems of general methodology for testing and data analysis were 
addressed using data from water tunnel tests. These tests were conducted 
with a 2.5% scale model of the F-16XL at Rolling Hills Research 
Company to identify nonlinear and unsteady aerodynamic model 
parameters. Although forced oscillation tests were performed about all 
three-body axes, only roll axis oscillation data were analyzed. 

Harmonic analysis of the small-amplitude oscillatory data was performed 
on single oscillations formed by averaging 30 cycles of data. This 
analysis provided a basis for the nonlinear analysis. Harmonic analysis 
was also applied to the large amplitude data. 

A general formulation of the aerodynamic model for aircraft that includes 
nonlinear unsteady aerodynamics was attempted. The resulting model 
could predict a broad class of aerodynamic responses of a rigid-body 
aircraft. The structure allowed easy interpretation of the model parameters 
by retaining conventional stability and control derivations for static and 
dynamic terms. Unsteady terms in the form of indicial functions were 
added to allow frequency dependent phenomena to be included. 

Experimental data and resulting model for the rolling-moment coefficient 
are presented in the ALAA Paper (see Publication 3). 

c) The oscillatory data obtained from the NASA LaRC 12-Foot Low Speed 
Tunnel for a 10% model of an F-16XL were used in previous years for 
mathematical modeling. 

A portion of the data; however, was analyzed for obtaining a measure of 
their accuracy. In this analysis the accuracy of the time stamp and the 
sample of time fidelity were first evaluated. Then, the effect of the timing 
signal was investigated. In the following step, the commanded and 
measured input signal was compared. The last part of the accuracy 
assessment was the evaluation of measured outputs. A measure of 
precision was obtained from repeated measurements under the same 
conditions. For the final statement in the accuracy, the outlier rejection 
rules were applied. The results of the data accuracy assessment are 
summarized in NASA TM Report (see Publication 4). 


2. Aerodynamic characteristics and instrumentation system calibration of Free 
Flving Aircraft Sub-scale Experimental Research Aircraft (LASER) 

The FASER program was established at the NASA LaRC to provide an 
affordable, easy-to-modify testbed to conduct research in the area of stability and 
control and system identification. Since the remote-control airplane will be used 
in high amplitude rapid maneuvers, it was desired to have a simulation program of 
FASER to allow the pilot to practice the maneuvers before the actual aircraft is 
flown. In addition, the simulation program can be used to design the experiment 
and to aid control law development. 

In this part of the research the existing aerodynamic model equations were 
updated using wind tunnel data. These data were obtained in the form of 
aerodynamic coefficients and their variation with the angle of attack, sideslip 
angle and control surface deflections. The analytical expressions for these 
dependencies were developed by applying the multivariate orthogonal modeling 
technique. An example of the resulting equations is given in Table 1. 

Wind tunnel data were also analyzed for obtaining calibrations of the alpha and 
beta-vanes. The analytical forms of these calibrations were developed by a 
stepwise regression. Examples of the calibration equations are given in Table 2. 

The complete results of the research conducted are in the NASA DSCB Internal 
Memo (see Publication 5). 

3) System Identification Methodology 

The past, present and future of system identification applied to aircraft at NASA 
LaRC were discussed in the Journal of Aircraft (see Publication 5). 

The purpose of this paper was to present a general approach to aircraft system 
identification, as formulated at NASA LaRC. The paper starts with a historical 
survey of flight testing and stability and control parameter estimation at NACA 
Langley Memorial Aeronautical Laboratory, and later at NASA LaRC, which lead 
to the present state of the art. The individual steps in the methodology are 
introduced and briefly explained. Four recent examples are mostly of the 
methodology application are given. These examples are mostly related to 
problems that go beyond traditional stability and control estimation from flight 
data. The paper is concluded by some thoughts about the future development of 
the methodology. 


Tablel. Aerodynamic model equations for two different subspaces 


-1.5° < a <10° and -10 < 0 <10 

C D (a,0,6 e ,8 f ) 

.0405 + 1.6365a 2 + .212a8 f + .2631a8 , - .041 5] 5 } + 
.062a + .01 3^ + .108 15} +.042 68 f 

C L (a,0,8 e ,5 f ) 

4.2977 a + .627 55 f +.291 58, + .046 -.6088} -2.%l\\cc l S f 


-32418, -.4917a -.0335/? + . 0646^ -.8846a 2 
+ .9092 0*8, +.2939; 6* -. 041106, -.1356oc8* f + 
.052 508 f + .06666* +.44556* ~.68690*8 { -.0294 

C y {a,0,8 a ,8,) 

-350 + .08936, - .02 56 a + .02965a - . 1 66206* - .04462 a8, +.11 a0 

C 1 (a,0,S a ,8,) 

-.198 8 a -.04110-. 221 a0 + . 245125] +.0212a8 a 

C n (cc,0,S a ,S r ) 

- .0435£ r + .0584/?-. 1 1 8 \x0+ 0164208} +.02352z<? r -.1350*8, 

- 34690? 0+.O1 266* 


10° < a < 20" and - 10 < /? < 10 

C D (a,0,8„8 f ) 

~.251 + 1.99a+.404a5, +.5586* + .1185} -.138} -.185 cc6 e 8 f 
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-11.892a 2 
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.348 -.8249^ -3.12ct+2.183a0* +132a* 8, +1.116a* +.8948^ 
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-310 + . 096198, +.81a8, -,417a 2 8 r + 3316a+ 381a0* ~.2O25a0 + 
3960*8, - .0898205*, - .00615] +15.41a 3 - .2164- .2620* ~.112l8 a 
-\.643a*8 a -12.7a 2 -.069 8} 

C l (a,0,8 a ,6 r ) 

-.131 8 a -9 .la* 0 - 368a8 a - .037 - .1591a0* -.008965} +1.1 31a* 5 a 
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-.004 08 a 8, -1.3773a 3 + .306a 
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-1.23a 2 0-maS} -.0007 8 a 8, -.888a 3 - .02836 -.083<? 
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Table 2. Calibration equations for alpha and beta vanes 


Region 

Relative Voltage 
Limits 

Equation 

-7.5° <a<10° 

or nv > .43681 

a" = 99.539 -201. 8 la„, - 4.2825 j3 nv 

10° <a<40° 

.29692 < a m < .43681 

a 0 = 106.98 - 229.71a,, + 45.279#, - 68.554#, 

40" < a < 80° 

a„ <.29692 

a° = 98.419 -180.03a„ + 40.879# -164.08a,,#, 

10° <^<30° 

Pm - -51086 

P° = -229.29 + 640.23#, - 35 1.85#, + 39.075a w 
-62.05a B 

^nvr'nv 

-10° </?<10° 

.42711 < (3 m < .51086 

P° = -11741.1 + 7222.2#, -15273#, +11031#, 
+ 60.448a,, - 170.78a 3 + 176.99a 3 , 


-30° < p <-10 


A, £.42711 


P° = -148.06 + 456.25 A, -326.72#, + 8.4967a„ 
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Figure 1. Effect of angle of attack and amplitude on lateral coefficient, C l , for rolling 

oscillations with k = 0.171. 



- 0 - k=0.073 
H- k=0. 1 22 
k=0.171 
k=0.19 
— - k=0.22 
-0- k=0.232 
-V- k=0.268 


a, deg 


Figure 2. Variation of in-phase and out-of-phase components of rolling moment 
Coefficient with angle of attack. Rolling oscillations. 
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Figure 3. Estimated time delay and gain parameter with their 2-sigma intervals 
Rolling-moment coefficient. Rolling oscillations. 
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Figure 4. Estimated rolling moment parameters from static and oscillatory tests 






PUBLICATIONS 


1. Klein, Vladislav; Murphy, Patrick C. and Szyba, Nathan M.: “Analysis of Wind Tunnel 
Lateral Oscillatory Data of the F-16XL Aircraft.” NASA TM-2004-2 13246. 

2. Szyba, Nathan M.: “Unsteady Aerodynamic Modeling of Dynamic Oscillatory Wind 
Tunnel Data of the 18% Scale F=16XL Aircraft.” M.S. Thesis, The George Washington 
University, 2004. 

3. Murphy, Patrick C.; Klein, Vladislav and Szyba, Nathan M.: “Progressive Aerodynamic 
Model Identification From Water Tunnel Test of the F-16XL Aircraft.” AIAA Paper 
2004-5277. 

4. Kim, Sungwan; Murphy, Patrick C. and Klein, Vladislav: “Evaluation and analysis of F- 
16XL Wind Tunnel Data from Dynamic Tests.” NASA TM-2004-2 13234, April 2004. 

5. Eissa, Mohamed M.: “Model Structure Definition and Analysis For the Wind Tunnel 
Force and Moment Coefficients For a Sub-Scale Research Airplane.” NASA DSCB 
Internal Memo, August 2004. 

6. Morelli, Eugene A. and Klein, Vladislav: “Application of System Identification of 
Aircraft at NASA Langley Research Center.” Journal of Aircraft, Vol. 41, No. 2, 
January-February, 2005, pp. 12-25. 


